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Abstract 

We study a new hypothetical form of solid carbon C28C, with a unit cell 
which is composed of the C28 fullerene cluster and an additional single car- 
bon atom arranged in the zincblende structure. Using ab initio calculations, 
we show that this new form of solid carbon has lower energy than hyper- 
diamond, the recently proposed form composed of C28 units in the diamond 
structure. To understand the bonding character of of these cluster-based 



solids, we analyze the electronic structure of C28C and of hyperdiamond 
and compare them to the electronic states of crystalline cubic diamond. 

PACS numbers: 61.46,+w, 61.66.Bi, 71.25.-s 

1 Introduction 

The C28 cluster has aroused considerable interest recently. This unit is the smallest 
fullerene that has been produced in significant quantities in experiments of laser 
vaporization of graphite |J. The structure of C 2 8 [0 is shown in Figure [l]. The 
twelve pentagons and four hexagons that comprise this structure are arranged in 
a pattern that gives rise to three nearly tetrahedral bond angles around four apex 
atoms (two such atoms are marked by A in Figure [I]). Dangling sp 3 orbitals on 
the four apex atoms render this cluster chemically reactive. A C28 cluster forms 
stable compounds when it is produced with U, Ti, Zr, or Hf atoms trapped at its 
interior gj, g, §. 

There are three inequivalent sites on the C28 cluster: we refer to atoms at 
the apex sites as cage-A atoms, their immediate neighbors as cage-B atoms, and 
the ones that are not connected by covalent bonds to the cage-A atoms as cage-C 
atoms (see Figure [I]). It has been proposed that the C28 unit could be stabilized by 
externally saturating the dangling bonds of the four cage-A atoms with hydrogen 
0) H §|> 01- Saturation of the dangling bonds through intercluster covalent bonding 
has also been considered: The four cage-A atoms of the C28 unit form the vertices 
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of a tetrahedron, making the C28 unit analogous to a tetravalent atom. A natural 
choice for a crystal composed of these units is the diamond lattice [jlj, |8|, [|. This 
hypothetical solid, called hyperdiamond and symbolized by (023)2, is shown in 
Figure |2|(a) in a perspective view. In Figure 0(b), the same structure is shown 
with all the cage-A, cage-B, and cage-C atoms on a single (110) plane indicated. 

A second candidate which allows for tetrahedral bonding of C28 units is a com- 
pound of C 2 8 clusters and tetravalent atoms in the zincblende structure [§]. Since 
the hyperdiamond lattice contains large voids in its structure, one may expect 
that a solid composed of a more compact packing of C 2 s clusters alternating with 
individual tetravalent atoms, would be energetically more favorable. This kind 
of structure can be obtained by simply replacing the central cluster in Figure ^ 
with a tetravalent atom and bringing the neighboring clusters closer to the cen- 
tral atom to form covalent bonds. In this paper we consider a carbon atom as 
the second component of the zincblende structure. We have performed extensive 
first-principles calculations to give a detailed comparison between this zincblende 
structure and hyperdiamond. The results are also compared to cubic diamond. 

The rest of this paper is organized as follows: Section describes the compu- 
tational approach used in our first principles calculations. Section |3| demonstrates 
the energetic favorability of C28C over (C2s)2 through total energy comparisons. 
In Section | the electronic structure of C28C is discussed in detail and compared 
to that of diamond and of (023)2- Finally, in Section [5] we draw conclusions on 
the new C28 based solids and comment on other possible solids based on similar 
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cluster-atom combinations. 



2 Methods 

We obtained total energies and electronic densities, as well as single particle elec- 
tronic states and eigenvalues by carrying out calculations within the framework of 
density functional theory and the local density approximation (DFT/LDA) JED), |TT|] . 
The optimal atomic coordinates and the equilibrium lattice constants are obtained 
with a planewave basis, including plane waves with kinetic energy up to 36 Ry. We 
also performed similar total energy calculations for bulk diamond, in order to make 
consistent total energy comparisons, and to obtain initial estimates for equilibrium 
bond lengths in C28C. For the reciprocal space integration, the sampling k-points 
have been chosen in a way such that the density of k-points in the first Brillouin 
zone is kept approximately constant for all the structures we considered. For bulk 
diamond, we used 125 k-points in the full Brillouin zone (which corresponds to 
eleven special points in the irreducible Brillouin zone) ||12|| , with correspondingly 



smaller sets for the C28C and (C2s)2 structures. Within the molecular dynamics 



framework of Car and Parrinello [13|, we used the steepest descent method for the 
initial relaxation of the electronic degrees of freedom, followed by the more efficient 
conjugate gradient method close the Born-Oppenheimer surface. 

The ionic potential, including the screening from core electrons, was modeled 
by a nonlocal norm- conserving pseudopotential |TJJ], and the Kleinman-Bylander 
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scheme |15| was employed to make the potential separable in Fourier space. The d 
angular momentum component was treated as the local part of the potential with 
the s and p components containing the nonlocal contributions. 

The choice of lattice constant for the C 28 C lattice was guided by the diamond 
calculations. In C28C we first chose the length of the bond between the single 
carbon atom and the C28 cluster to be the equilibrium bond length of diamond, 
as obtained by our calculations. Since the interaction between cage-B atoms of 
neighboring C 2 8 clusters can affect the optimal length of this C-C 2 8 bond, we 
performed the calculations at several different lattice constants in the range near 
our initial choice, corresponding to different values for the C-C28 bond length, 
while the internal structure of the C28 cage was held fixed at the one determined 
by hydrogenation of the four cage-A atoms of the isolated cluster. We used this 
structure because the hydrogen atoms saturate the dangling bonds of C 28 just as 
the single carbon atoms do in C28C. By fitting to a Birch-Murnaghan equation of 



state [IS], we found that the lowest energy lattice constant corresponds to a C-C28 
bond length of 1.53 A. Fortuitously this was one of the actual lattice constants for 
which the C 2 sC calculations were performed. For the optimal lattice so obtained, 
the positions of the atoms were then relaxed by minimizing the magnitude of 
forces obtained through the Hellmann-Feynman theorem. The cutoff below which 
the forces were considered to be negligible is 0.01 Ry/a.u. For the case of C 2 sC, 
this relaxation not only lowers the energy per atom but also changes the electronic 
structure from a metal to a semiconductor, which makes the crystal more stable. 



We discuss this further in Section |j. Similar calculation procedures were performed 
to obtain the optimal structure in (C2s)2- 

3 The Relative Total Energies 

The optimal total energy per atom for C28C and (C2s)2 after full relaxation are 
0.453 eV and 0.744 eV respectively, where we have used the energy per atom of 
diamond as the reference. 

The energy of C28C is lower than that of hyperdiamond, which is due in part to 
the presence of the additional fourfold coordinated carbon atoms (the added single 
carbon atoms). However, intercluster interactions also play an important role in 
stabilizing C28C. We can get an estimate of the effects of intercluster interaction 
by simply comparing the energy between a unit cell of C28C and the equivalent of 
twenty-eight atoms of (C2s)2 plus one atom of diamond: the C28C lattice is lower 
in energy by 0.26 eV per atom. This should be compared to the energy difference 
between C28C and (C2s)2 quoted above, which is 0.291 eV. Thus, intercluster 
interactions are more favorable in C28C, rendering it lower in energy than it would 
be if the single added atoms were equivalent to diamond atoms and the C 2 8 clusters 
were equivalent to those of (C 2 8)2- 

We discuss next the effect of C28 cage relaxation on our results. The relax- 
ation has decreased the energy per atom for the case of C28C by as much as 0.24 
eV/atom, while only approximately 0.01 eV/atom is obtained through relaxation 
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from the same original choice of cage (that of the hydrogenated C28 cluster) in 
the case of (C2s)2- In order to understand this, we first compare the distances 
between atoms in neighboring clusters between different structures and how re- 
laxation affects them. In hyperdiamond, other than the covalent bond between 
cage-A atoms on neighboring clusters, the smallest distance between two atoms on 
neighboring clusters is between a cage-A atom on one cluster and a cage-B atom 
on the next cluster. This distance is 2.57 A. The next smallest distance between 
atoms on neighboring clusters is between two cage-B atoms and has a value of 3.15 
A [see Figure 2(b)]. In this case, relaxation does not induce any significant change 
on the geometry of the cage structure relative to the structure of the isolated, 
hydrogenated C28 cluster. 

In contrast, in C 2 sC the closest distance between atoms on neighboring clusters 
before relaxation is only 2.00 A, between two cage-B atoms of the original cages. 
After the relaxation, this distance is reduced to 1.64 A, which evidently introduces 
additional bonding between cage-B atoms as we will see from the analysis of the 
electronic states in Section Therefore, as far as the structural relaxation is 
concerned, we have observed significantly different behavior between C28C and 
(C28)2- This was expected because the closest intercluster distance in the C28C 
solid occurs between two cage-B atoms, both of which are threefold coordinated. 
An additional B-B intercluster bond will be energetically favorable as carbon atoms 
prefer to be fourfold coordinated in this environment. In the case of (C2s)2 the 
closest intercluster distance is between a cage-B and a fourfold coordinated cage-A 
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atom, which makes the original cage structure in (C2s)2 electronically more stable 
compared to C28C. 

In Table [I] we display the bond lengths of C28C and (C2s)2 before relaxation (the 
hydrogenated structure) and after relaxation. The relaxation in (C 2 s)2 produces 
lengthening of the intercluster bonds between cage-A atoms on neighboring C28 
units and shrinking of all other intracluster bonds. The intracluster bonds which 
are shortened by the largest amount are those between cage-C atoms, which are 
the farthest bonds from the cage-A atoms. In the case of C28C, a somewhat 
different relaxation pattern of the cage geometry emerged: the C-C28 bond and 
the C-C intracluster bonds were shortened while the other two intracluster bonds 
were lengthened. This relaxation can be attributed to the formation of additional 
inter-cluster bonds between cage-B atoms. 

4 Electronic States 

In Figure (a), (b), and (c), we display the total valence electron densities of 
the fully relaxed C28C and (C2s)2 solids and diamond on the (110) plane. The 
length scales in Figure ^| (a) and Figure ^| (b) have been chosen so that these plots 
cover approximately equal areas. It is evident from a comparison of these figures 
that C28C has a higher atomic density than (023)2- Is is also apparent that the 
basic C28 cluster geometry remains essentially unchanged in both the C28C and 
(028)2 structures. This observation assures us that the C28 units are not altered 
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significantly through the introduction of other C28 clusters or additional carbon 
atoms in the solid forms. The charge density around the additional carbon atom 
outside the cage (denoted by the symbol X in Figure 0(a)) is seen to be similar 
to the intercluster A-A bond in (C 28 ) 2 . These bonds are single covalent bonds 
between two fourfold coordinated carbon atoms: comparison to Figure 0(c) 
shows that they have the same bonding character as the bonds in bulk diamond. 
The main difference between (C 28 ) 2 and C 28 C is the presence of additional bonds 
between cage-B atoms of neighboring clusters in the C 28 C crystal. These bonds, 
seen clearly in Figure 0(a) between the two cage-B atoms of neighboring clusters, 
are somewhat weaker than regular covalent bonds, as expected from the fact that 
they are 1.64 A long, compared to the 1.54 A bond distance in bulk diamond. 

In Figure f| we display the density of states (DOS) of C 28 C and (C 28 ) 2 at the 
optimal lattice constants with full relaxation. The DOS of diamond, as calculated 
in this work, is also shown in dashed lines for comparison. The DOS of the cluster- 
based solids exhibit many features that can be related to features of the diamond 
DOS. For instance, the total valence band width in all three cases is essentially 
the same, 21 eV. Some specific features deserve closer attention: In both cluster- 
based solids, the states corresponding to the intercluster bonding between cage-A 
atoms, which is similar in nature to the bonding in diamond, are far below the 
fermi level. States immediately below and above the fermi level derive from the 
bonding properties of atoms within the clusters. A careful examination of the 
symmetry of the wave functions of the highest occupied and lowest unoccupied 
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states indicates that the states of C28C immediately below the fermi level are due to 
a combination of intercluster bonding states between cage-B atoms and 7r bonding 
states between cage-C atoms within the cluster. The same type of state involving 
intercluster bonding among cage-B atoms is responsible for the highest occupied 
state in (028)2, even though the interaction is much weaker in that case. On the 
other hand, the lowest unoccupied states in both solids are primarily due to a 
combination of intercluster bonding states between cage-B atoms and antibonding 
states between cage-C atoms within the cluster. 

The DOS reveals that (C 2 s)2 is a semiconductor with direct band gap, which 
is consistent with results reported for (C2s)2 in previous works |], H fl- The 
fully relaxed C28C structure is also a semiconductor with direct band gap equal 
to 1.16 eV, approximately two thirds of the band gap of (C 2 8)2, and much smaller 
than the band gap of diamond. As we have noted above, there is a much stronger 
intercluster interaction in the case of C28C than in (C2s)2; it is therefore reasonable 
to expect that the electrons should be more delocalized in C28C than in (C2s)2, 
which leads to the smaller band gap of C28C. Here we wish to remind the reader 
of the well known inability of DFT/LDA to reproduce the experimental band gaps 



in semiconductors and insulators [|17]], so that all the numbers quoted above are 
underestimates (by approximately a factor of 2) of the true band gaps, if these 
solids would be realized. For example, the band gap of crystalline cubic diamond 
obtained by the present calculation is 3.85 eV, which is approximately 2/3 of the 
experimentally measured value 5.52 eV. 
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5 Conclusions 

In summary, we have studied two cluster-based solid structures, the C28C and 
(C 28 ) 2 lattices. Both structures satisfy the condition that the dangling orbitals 
with sp 3 character on the cage-A atoms are fully saturated. Using first principles 
calculations, we demonstrated that C 28 C is energetically preferred to (C 28 ) 2 . The 
crucial role of intercluster interactions at the cage-B atoms is revealed by the 
structural relaxation in C 28 C, energetic comparisons, and the valence electronic 
charge density. The electronic structures of both C 28 C and (C 28 ) 2 were analyzed. 
It was shown that for both solids the intercluster bonding between cage-B atoms 
and 7r bonding between cage-C atoms within the cluster contribute to the highest 
occupied states in C 28 C. The lowest unoccupied states, on the other hand, come 
from a combination of the cage-B type intercluster bonding states mentioned above 
and anti-bonding states between cage-C atoms within the cluster. 

Just as intercluster interactions render C 28 C lower in energy than (C 28 ) 2 , other 
solids based on a C 28 unit with saturated dangling bonds might be candidates for 
stable structures as well. For example, we expect that solids analogous to C 28 C, 
but with the single C atom that links the C 28 units replaced by other group-IV 
atoms (Si, Ge, Sn) would be equally stable, and with similar electronic properties. 
Alternatively, solid forms composed of C 28 units connected by O atoms (chemical 
form C 28 2 ) in a manner analogous to silica or crystalline forms of Si0 2 could 
be rather stable and somewhat easier to make, given the flexibility of packing of 
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tetrahedra bonded at their corners. The conditions under which such solids may 
be experimentally realized remain to be investigated. 
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Table 1: Bond lengths before and after relaxation of the cage structure for both 
the (C 2 8)2 and C 2 sC structures. The three types of intracluster bonds (bond A-B, 
B-C and C-C) refer to the labels of atoms shown in Figure 1 and Figure 2(b). 

Figure 1: Structure of the C 2 8 unit. The dashed lines indicate two of the four 
axes of C3 symmetry; the dashed-dotted line indicates one of the three axes of C2 
symmetry. The three atom types are: (i) an apex or cage-A atom , where three 
5-fold rings meet, (ii) a cage-B atom which is part of a 6-fold ring and is bonded 
to a cage-A atom, and (iii) a cage-C atom which is also part of a 6-fold ring, but 
is not bonded to a cage-A atom. 
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Figure 2: (a) Perspective view of the structure of the hyperdiamond lattice. The 
C 28 clusters are bonded together at the cage-A atoms shown here in black, (b) The 
same structure as (a) shown along a (110) crystallographic direction, where the 
positions of cage-A, cage-B, and cage-C atoms are indicated on a single plane. The 
dashed line indicates the place where weak intercluster B-B type bonding occurs in 
(C 28 ) 2 . The C 28 C structure is obtained by replacing the central cluster by a single 
C atom and bringing the neighboring clusters closer to that site to form covalent 
bonds. In that case the interaction between neighboring cage-B atoms becomes 
much stronger (see text). 

Figure 3: Total valence electron density on the (110) plane of (a) C 28 C, (b) (C 28 ) 2 
and (c) diamond. The white areas represent the highest electron density and the 
black areas the lowest. The positions of cage-A, cage-B and cage-C atoms are 
shown in (a) and (b). The positions of the single extra carbon in C 28 C is indicated 
by the symbol X; the same symbol indicates the position of atoms in the diamond 
lattice (c). 

Figure 4: Density of states vs. energy for the C 28 C (upper part) and (C 28 ) 2 (lower 
part) structures. For comparison the density of states of diamond (as calculated 
here) is also shown in dashed lines. The fermi level (top of the valence band) has 
been taken to be the zero of the energy scale for each structure. 
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